Slow modulations of periodic waves in Hamiltonian 
PDEs, with application to capillary fluids 

S. Benzoni-Gavage* P. Noble} and L.M. Rodrigues* 
Universite de Lyon, CNRS UMR 5208, Universite Lyon 1, 
Institut Camille Jordan, 43 bd 11 novembre 1918 
F-69622 Villeurbanne cedex 

March 27, 2013 



Abstract 

Since its elaboration by Whitham, almost fifty years ago, modulation theory 
has been known to be closely related to the stability of periodic traveling waves. 
However, it is only recently that this relationship has been elucidated, and that 
fully nonlinear results have been obtained. These only concern dissipative systems 
though: reaction-diffusion systems were first considered by Doelman, Sandstede, 
Scheel, and Schneider [Mem. Amer. Math. Soc. 2009], and viscous systems of con- 
servation laws have been addressed by Johnson, Noble, Rodrigues, and Zumbrun 
[preprint 2012]. Here, only nondissipative models are considered, and a most ba- 
sic question is investigated, namely the expected link between the hyperbolicity of 
modulated equations and the spectral stability of periodic traveling waves to side- 
band perturbations. This is done first in an abstract Hamiltonian framework, which 
encompasses a number of dispersive models, in particular the well-known (general- 
ized) Korteweg-de Vries equation, and the less known Euler-Korteweg system, in 
both Eulerian coordinates and Lagrangian coordinates. The latter is itself an ab- 
stract framework for several models arising in water waves theory, superfluidity, and 
quantum hydrodynamics. As regards its application to compressible capillary fluids, 
attention is paid here to untangle the interplay between traveling waves/modulation 
equations in Eulerian coordinates and those in Lagrangian coordinates. In the most 
general setting, it is proved that the hyperbolicity of modulated equations is indeed 
necessary for the spectral stability of periodic traveling waves. This extends earlier 
results by Serre [Comm. Partial Differential Equations 2005], Oh and Zumbrun 
[Arch. Ration. Mech. Anal. 2003], and Johnson, Zumbrun and Bronski [Phys. D 
2010]. In addition, reduced necessary conditions are obtained in the small amplitude 
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limit. Then numerical investigations are carried out for the modulated equations of 
the Euler-Korteweg system with two types of 'pressure' laws, namely the quadratic 
law of shallow water equations, and the nonmonotone van der Waals pressure law. 
Both the evolutionarity and the hyperbolicity of the modulated equations are tested, 
and regions of modulational instability are thus exhibited. 

Keywords: Whitham modulated equations, traveling wave, spectral stability, modulational 
instability, Lagrangian coordinates. 

AMS Subject Classifications: 35B10; 35B35; 35Q35; 35Q51; 35Q53; 37K05; 37K45. 



1 Introduction 

This work is motivated by the nonlinear waves analysis of the so-called Euler-Korteweg 
system, which arises in the modelling of capillary fluids - these comprise liquid-vapor 
mixtures (for instance highly pressurized and hot water in nuclear reactors cooling system, 
in which the presence of vapor is actually dramatic), superfluids (Helium near absolute 
zero), or even regular fluids at sufficiently small scales (think of ripples on shallow water 
or other thin films). In one space dimension, the most general form of the Euler-Korteweg 
system we consider is 



(1) 



d t p + d x {pu) = 0, 

d t u + ud x u + d x {E p S) = 0. 



in Eulerian coordinates, its counterpart in mass Lagrangian coordinates being 
(2) 




8yU, 

d y (E v e) , 

where p is the fluid density, v — 1/p its specific volume, u its velocity, viewed either as 
a function of (t,x) or as a function of (t, y), where y is the mass Lagrangian coordinate 
(by definition, dy = pdx — pudt). The energy density $ in ([I]) and the specific energy e 
in ([2]) are related through § = pe, or equivalently e = vd?, and are regarded as functions 
of (p,Px) and (v,v y ) respectively. In those systems, the notation E stands for the Euler 
operator, that is 

d<S „ ( dS\ _ d<e _ (Be 



E P^ = ~x D z a~ ' = o D y 

Op \VPxJ <JV \UV y 

where D x and D y mean total derivatives. A widely used class of energies, dating back to 
Korteweg's theory of capillarity, read 

(3) *(p,p*) = f(p) + \ jr(p)(p x )\ 
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or equivalently 

(4) 4y,v v ) = f{v) + -K(v)(v y ) 2 , 

with the relationships F = pf, k = p 5 J^, and with various choices of (as far as is 
concerned, k is often chosen to be constant, whereas quantum hydrodynamics equations 
correspond to p,Jfc = constant). 



Abstract framework Equations ([TJ and ^ fall into the class of abstract Hamiltonian 
systems of evolution PDEs of the form 

(5) d t u = a ?(Eje>\y\), 

where J? = d x J is a skew-symmetric differential operator, J being a symmetric, nonsin- 
gular matrix with constant coefficients, Jf? is a functional involving first order derivatives 
only, and E denotes again the Euler operator: 

E^[U] a = ^(U,U x ) - D x (J^(\J,U x ^j , ae{l,...,N} 

if U has N components. Equation (|5| being space-invariant, it admits a conserved quan- 
tity called an impulse, say such that 

and since J is nonsingular we can explicitly take for £} the quadratic quantity 

«g(U) := |U- J _1 U. 
For further use, let us mention that associated with is the local conservation law 

(6) dt£(U) = d x {y\\j]) 

satisfied along any smooth solution of (|5]), where 

(7) y[u] ■.= u • E^r[u] + i_^[u] , ljt[u] : = u a , x ^— (u, u„) - jf(u, u x .) . 

The dot • here above is just for the 'canonical' inner product U • V = U a V a in R N , and 
the letter L stands for the Legendre transform (even though it is considered in the original 
variables (U, TJ X )). In its definition we have used Einstein's convention of summation over 
repeated indices, and we shall do so repeatedly in the sequel. 
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Examples The Euler-Korteweg system (JI]) fits into this framework with 

U=f£J, M>= ^pu 2 + g(p, Px ), £ = pu, y = p£ p <g + L£ ) 
as well as the system in Lagrangian coordinates (|2]), with 

U=^J, 34? = \i 2 + 4y,v y ) , £! = vu, S fi = vE v £+Le. 

An even simpler example is the scalar equation 

(8) d t v = d x (E v Ji?) , 

with j£? = |t> 2 , comprising the generalized Korteweg-de Vries equation (gKdV), 

(9) d t v + d x p(v) = -8%v , 



if we take 



j?=m + \vi, f = - P . 



In this case, 5? = —f — pv — vv xx + ^v 2 - and (|6| is the well-known conservation law for 
v 2 when v is solution to (gKdV). 

We will focus on the 'generic' situation in which the abstract system ^ admits a 
family of periodic traveling wave solutions parametrized, up to translations, by their 
speed and N + l constants of integration. This rather large number of degrees of freedom 
makes the stability of periodic waves a difficult problem. Furthermore, in the context of 
Hamiltonian PDEs we can only hope for neutral stability, namely that the full spectrum 
of linearized equations lies on the imaginary axis. On the one hand, this kind of spectral 
stability makes the nonlinear stability analysis more delicate to tackle than in situations 
involving dissipation processes that are likely to push the spectrum - except for the null 
eigenvalue linked to translation invariance - into the left half plane. (Nonlinear results 
have recently been obtained concerning reaction-diffusion systems [7], and viscous systems 
of conservation laws [IS]-) On the other hand, the underlying variational framework can 
be of great help to prove nonlinear stability results. This has been done for a wide 
variety of solitary waves, see [1] and references therein. The literature on the stability of 
periodic waves is much more limited, and nonlinear results are limited to stability under 
perturbations of the same period p]. Up to our knowledge, spectral stability has been 
proved for periodic wave solutions to the (standard) Korteweg-de Vries equation by using 
its integrability to compute explicitly the spectrum [1] , and for small amplitude periodic 
wave solutions to a limited number of dispersive equations, comprising the nonlinear 
Schrodinger equation and the generalized Korteweg-de Vries equation [10] . It is still a 
wide open problem for large amplitude periodic wave solutions to more general, dispersive 
PDEs. 

As regards spectral stability, it is more tractable if sideband perturbations only are 
considered. There are indeed several analytical tools - like for instance Evans functions 
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and associated winding numbers - to look for possible spectrum in the vicinity of zero. 
A related topic is modulational stability, in connection with Whitham's approach of slow 
modulations to periodic waves (the most famous modulational instability being the one 
shown by Benjamin and Feir for Stokes water waves). 

By - at least formal - asymptotic analysis, we can see that slow modulations to 
periodic waves are governed by a system of averaged equations. Its size is the dimension 
of the periodic orbits' manifold, here N+2. This will be made more precise in Section[2j of 
which the main purpose is to extend to our abstract framework a result previously known 
for (gKdV) jH] and for viscous systems of conservation laws [T71 [18], [19] . This result gives 
a quantitative relationship between the sideband stability of periodic traveling waves and 
the spectral properties of modulated equations. In particular, it shows that a necessary 
condition for the stability of a given periodic traveling wave is the (weak) hyperbolicity 
of modulated equations at the corresponding point in parameters' space. In other words, 
we give a rigorous proof that modulational stability is necessary for spectral stability, a 
result that is often taken for granted in the physics literature - back to 1970, Whitham 
himself was indeed saying that 'the relation of the stability of the periodic wave with the 
type of the [modulated equations is] given in the previous papers' [20]. In addition, we 
investigate in some detail how the modulated equations degenerate in the small amplitude 
limit, and receive a reduced system coupled with a 2 x 2 system for the wavenumber and 
the amplitude of the wave. This extends to our abstract framework observations that 
were made by Whitham [21]. Unsurprisingly, when applied for instance to the Euler- 
Korteweg system, that limit gives as a reduced system the lower-order, Euler equations. 
Hence a stability condition for small amplitude periodic waves: the Euler system must be 
hyperbolic at the mean value of the wave. This rather natural condition does not seem 
to have been pointed out earlier. A reason is certainly that modulated equations have 
mostly been considered for scalar models, like KdV or the Klein-Gordon equation in |21j . 
for which the hyperbolicity of the reduced model - the inviscid Burgers equation for KdV, 
the wave equation for Klein-Gordon - is trivial. 

In Section [3j we concentrate on the Euler-Korteweg system. We derive modulated 
equations in both kinds of coordinate systems, namely the Eulerian one ([TJ, and the mass 
Lagrangian one (|2]). In addition, we point out a nice - if not surprising — relationship 
between them. To be precise we show that, away from vacuum, the modulated system 
for ^ is equivalent to the modulated system for Q through a mass Lagrangian change 
of coordinates, hence the following commutative diagram. 





mass Lagrangian 






change of coordinates 




(!) 




(2) 


Whitham's . 
averaging 




I 




(1) 




(2) 



Then in Section [4] we go further into specialized cases and investigate in more details 
the periodic orbits' manifold when the energy is of Korteweg type (|3]). In this case, the 



5 



nature of the phase portrait associated with profile equations for traveling waves highly 
depends on the monotonicity and convexity properties of pressure in terms of volume [3]. 
We have considered two types of pressure laws, each one being motivated by a specific 
physical application. At first, we have taken a quadratic pressure with respect to density, 
which corresponds to shallow water equations - the Korteweg part of the energy then 
taking into account surface tension on the water surface. A more involved case arises 
with van der Waals type pressure laws, typically corresponding to liquid-vapor mixtures 
(with capillarity effects). In both cases - shallow water and van der Waals - , we have 
investigated numerically the hyperbolicity of Whitham's equations associated with various 
families of periodic waves. We have found rather large regions of hyperbolicity. Failure of 
hyperbolicity occurs for waves 'close to' unstable constant states/solitary waves or/and 
with sufficiently large periods. 



2 Periodic waves and modulated equations 
2.1 General material 

Traveling wave solutions A traveling wave solution to ^ of speed c is characterized 
by the profile equations 

(10) E(J4? + cJ?)[U] = A , 

where the components X a of A e Mr are merely constants of integration. As was observed 
by Benjamin [2], the impulse flux 5? is involved in the Hamiltonian associated with the 
Euler-Lagrange equations (10), which turns out to be 5? + c£l. In other words, solutions 
of (10) must satisfy 

(11) y[v} + c«s(u) = n 

for some new (scalar) parameter fi. The profile U may be viewed as a stationary solution 
of the abstract Hamiltonian system ^ rewritten in a frame moving with speed c, that is 

(12) d t \J = f{E{M> + c£)[\J}). 

We may also note that in this moving frame the additional conservation law ^ reads 

d t £(u) = d x ({y + c£)[\j}), 



which obviously admits the profile U as a stationary solution, according to (11). 



Linearized problem In order to investigate the stability of U as a solution of (12), we 
start by linearizing this system about U, which yields 

(13) d t \J = jzAJ , st := J? A , A := Hess(je + c£)[\l] . 
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The Hessians here above are given by 



Hess^[U] = J" 1 

whatever U, and 



d 2 3f d 2 3f ( d 2 Jf? d 2 

where all second order derivatives of J%? are evaluated at (U, TJ X ). By differentiating the 



profile equations (10) with respect to x we observe as usual with translation-invariant 
problems that AU^, = 0, hence also £^]J_ X = 0. Furthermore, we have several parameters, 
namely the speed c, and the constants of integration X a , /i. If they are all independent, 
as is typically the case with periodic traveling waves, we receive remarkable identities by 



differentiating the profile equations (10) with respect to those parameters. This yields 



X 1 



AU C = -E^[U] , hence &TU C = -U 
AU^ = , hence sfU^ = , 

AU Aa = e a , hence ^/U Aq = . 

(Here above we have denoted by e a the a-th vector of the 'canonical' basis in the U-space, 
and the subscripts c, fi, X a stand for partial derivatives with respect to those parameters.) 
When U is periodic, say of period 5, this period of course depends on the parameters c, 
X a , [A, and thus the derivatives U c , U , U Aq have no reason to be periodic. However, as 
pointed out in [6j, we can set up the generalized kernel of srf in the space of S-periodic 
functions with linear combinations of those derivatives. As a matter of fact, U^. and 
$ a := H^Hu — S (( U Aa are all in the S-periodic kernel of while 

* Q := {2, (U)} Aa , c U M + {S, (U)}c,,U Xa + {E, (U>W«U C 

is also S-periodic, and such that 

^* a = -{S,<U» A1 , Aa U ai eker^. 

Here above we have used the same convenient notation as in [6] 

{/) g}a,b = fadg ~ fb9a , 

and the brackets (■) stand for mean values on a period. The periodicity of the linear 
combinations mentioned above follows from the identity 

U a (H)-U Q (0) =S a U x (0), 

which holds true whatever the parameter a. Note that we find in this way at most N + 2 
independent elements of the generalized kernel of srf in the space of S-periodic functions 
if N is the dimension of the U-space, even though we have used iV + 3 candidates. 
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Modulated equations Following Whitham's 'two-timing method' [20], we search for 
solutions of ^ having an asymptotic expansion of the form 

U(t,x) = \J°(et,ex, (j)(et,ex)/e) + s\J 1 (et,ex, 4>{et,ex)/e,e) + o(e) , 

where U° and U 1 are 1-periodic in their third variable (subsequently denoted by 6). 
Denoting by X = ex, T = et the rescaled space variable and time respectively, 

k := (f>x , oj := 4>t , c :— —u/k , 

the very existence of a twice differentiable phase requires that 

(14) d T k + d x (ck) = 0. 

Plugging the asymptotic expansion in ([5| and using that d t = edx + uido, d x = edx + kde, 
we formally receive the following equations 



(15) 



d e (G° + cV°) = 



(16) <9 T U° = Jd x G° + JkdeiG 1 + cV 1 ) 

where the components G^ 1 , V®' 1 of G ' 1 and V 0,1 are given by 



V 



dig 



or equivalently V c 



j- x u° 



d 2 je TTl d 2 3f jaTTl ,^ / <9 2 ^f rrl <9 2 ^r , orrl 



dU a dU 



3U n dU, 



aUU/3 

,x 



dU a>x dUp~P ' dU a>x dUp jX 



dU n dU, 



,x 



dU a ~ 



kD, 



dUa^dUfc 



dxUj 



d 2 £ 



Vl := — — Ul , or equivalently V 1 = J^U 1 

OU Q Up 



all derivatives of being evaluated at (U°, kdeXJ ), and those of =S at U°. The zeroth 



order equation in (15) yields (G° + cV°) = constant, which just amounts to the traveling 



wave profile equation (10). More precisely, (G° + c V°) = A requires that, at fixed (T, X), 



U(x) := U°(T, X, kx) solves (10). The first order equation in (16) yields the averaged 
equation over 9 G [0, 1] 



(17) 



d T (U°) = 3d x (G°). 
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Together with the compatibility equation (14), this is the main set of modulated equations. 



There is an additional one associated with translation invariance of the original system. 
The fastest way to obtain it is to use the conservation law in ([6]), which yields 



(18) 
where 



d T (Q°)=d x (S ) 



Q°;=lV°-U , S° 



Xjo . G o _ ^ (U ; k9gV 0) + kdgU o ; 

OU a x 



the derivative of Jif being again evaluated at (U°, kdeU ). The other way is the main 



reason why we have written down G 1 in details. As a matter of fact, Equation (18) can be 
obtained by averaging the inner product of (16) with V°. The first term in this operation 
is simply 

(V° • d T U°) = d T (Q°) . 



On the right-hand side of (V° • (18)), we first have 



(V° • 3d x G°) = (U° • d x G°) . 
In order to deal with the other term, we can split G 1 + cV 1 into 

G 1 + cV 1 = A{kd e ) U 1 + C{kd e ) d x U° 
where A(kdg) and C(kde) are both linear differential operators with coefficients that are 



functions of (U°, kdgXJ ), and moreover Equation (15) means that 

A{kd ) dg\J° = . 

This implies that 

(V° • Jkde{A{kd e ) U 1 )) = - (U 1 • A(kd e ) kd e U°) = , 

hence 

(V° ■ JkdeiG 1 + cV 1 )) = - (kd e V° ■ C(kd e ) d x V } . 
Therefore, it remains to check that 

(U°-d x G°) - (kd e U o -C(kd )d x U°) = d x (S°), 

or equivalently, 

(G° • d x XJ°) + d x (kd e U° ■ V Vx J^ - 3>e) + (kd e XJ° ■ C(kde) d x XJ°) = . 
Now, recalling the definition of G° and making an integration by parts we get 
(G°-^U°) = ((a x U°).Vu^+(^U ).V u „jr) = d x (J?) - ((d^devyv^) 
= dx(JF - kd e U° ■ Vu x J4f) + kd x ((d e U°) ■ V Vx J^) . 
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Finally, recalling the definition of C and making an integration by parts, we have 



(deU -C(kde)d x U°) = 

-wish -<***(!£)> 

thanks to another integration by parts. 



Low frequency analysis We assume that U is a given periodic traveling wave profile 
of period S = 1/k, and that the set of nearby periodic traveling wave profiles is a iV + 2 
dimensional manifold if N the dimension of the U-space. As explained above, natural 
parameters for this manifold are the speed c of waves, the constants of integration X a 



for a G {1, . . . , N} as well as // showing up in the profile equations (10)-(11). In fact, 
we shall prefer a parametrization that is more natural in connection with Whitham's 
modulated equations, and assume that the manifold of periodic traveling wave profiles is 
parametrized by their wave number k (inverse of period), and the mean values 

(19) M=(U), P=<J2(U)). 



For a discussion of this assumption, see Appendix |B.2| In addition, we rescale all periodic 
profiles into 1-periodic ones, in such a way that k appears explicitly in the profile equations. 
For simplicity we still denote by U the profile now viewed as a function of 9 = kx, and 
similarly any nearby profile U is viewed as a function of 9 = kx. The latter must therefore 
satisfy 

(20) do — (U, kd e V) - kdo ^(U, kdoV) + c — (U) = , a G {1, . . . , N} . 



\dU a ' \dU a , x ' J dU a 

Our previous assumption means that for all (k, M, P) close to (k, M, P) there is a unique 
speed c = c(k, M, P) and a unique profile up to translations U = U(6 ) ; k, M, P) that is 
1-periodic, closd^to U, and solution to (19)-(20). For simplicity again, we just denote by 



A the differential operator A(kdg) considered in the previous paragraph with U° replaced 
by U (which amounts to Hess(^f + ci£?)[U] where d x is replaced by kd$), and by srf the 
operator 

£?{kd e ) = kd e JA(kd e ) . 

More explicitly, we have 

8 2 W d 2 W f d 2 W d 2 W 

dU a dU/3 dU a dUp )X \oU^ x dU^ oU a>x dU^ x 



1 We say that two 1-periodic functions are close to each other if their distance with respect to the sup 
norm up to translations is small. 
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Equation (20) equivalently reads AdeU = 0. Besides, this equation implies that 



d e (AU a + c a 3 1 U) =0, hence £/U a = -c a kd e U, 

where the subscript a stands for a partial derivative with respect to any parameter among 
M a and P. According to our assumption on the parametrization of periodic profiles, this 
makes at least N + I independent elements of the generalized kernel of s$ in the space 
of 1-periodic functions, and in fact N + 2 counting dgXJ. Let us mention straightforward 
identities for the formal adjoint 

£/* = -AJkd e . 

We clearly have indeed 

(Recall that e a is just a constant vector in the U-space.) Now, we are not only interested 
in the spectrum of srf in the space of 1-periodic functions but in the whole space of 
bounded functions. This is why we introduce the Bloch operators 

st v ■= &?(k(d e + iv)) , v G M/2ttZ , 

of which the spectra in the space of 1-periodic functions give the one of s$ on L°°: 

<j{st) = [J a{sf). 

i/eM/27rZ 



2.2 Modulational stability vs spectral stability 

Given the material introduced in §2.1[ we can show the following. 



Theorem 1. Let us assume that U is the profile of a periodic traveling wave solution to (|5]) 
of period 1/k and speed c, and that the set of nearby periodic traveling wave profiles U of 
speed c close to c is a N+2 dimensional manifold parametrized by (k, M = (U), P = (Q)), 
where 1/k is their period and 

g^iu-J^u. 

Let us consider the modulated system 

d T k + d x (ck) = 0, 
(21) { d T (U) = 3d x (G), 



where 



d T (Q) = d x (S) 



G a := ^f(U,kd e V) - kd e (^(\J,kd e U) 

dU a \oU a , x 

S := U • G - JT(U, kd e U) + (kd e U a ) kd e U) . 
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We also assume that the generalized kernel of srf in the space of 1-periodic functions is of 
dimension N + 2. Then a necessary condition for U to be stable is that the system (21) 
be 'weakly hyperbolic' at (k, M, P), in the sense that all its characteristic speeds must be 
real. 



Proof. It is based on a perturbation calculation, which relates the matrix of (21) at 
(k, M, P) to the one of srf u restricted by spectral projection to a (N + 2) dimensional 
invariant subspace. We first introduce the expansion 

where = £/° = JkdgAS ' is just srf viewed as an operator acting on 1-periodic 

functions, as well as A^ ^ is just A acting on 1-periodic functions, 

fit® := J(A(°) + kd e A^) , ^ : = J(A« + kd e A^) , ^ := JA& , 

dU a dUp tX dU^dUp \dU a)X dU P)X J dU a<x dU^ x 

(A^V) ■= Q2 '^ V, 



Differentiating (20) with respect to k we find that 

(22) ^/ (0) U fc + ^ w d e V = - c k kd 6 V . 

Indeed, we find at once that for all a G {1, ... , N}, 

3, (,A<». Ut + Ck J-U). + ^ m - D, (^) - tD. S .U t 
hence 

)U fc + c fc kd e V + 3kd e A {1) d e V = , 



which is equivalent to (22) since AdgXJ = 0. 
As already mentioned, 

$° :=fcU, $° :=U Ma , *; +1 :=U P 

all belong to - and span - the generalized kernel of while 

*° :=e aj ae{l,...,JV}, < +1 := J^U 

belong to - and span - the kernel of srf* . By definition of the mean values M a and P in 
(19), those functions are such that 

« ■ *?> = <w • 
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Therefore, we can add in a function \I>q such that (\I/q, vE^, . . . , ^° N , ^° N+ i) be dual to the 
basis ($Q, ■ ■ ■ ) ^°n+i) °f the generalized kernel of and span the generalized 



kernel of stf* . Recall that by our main assumption, is an isolated eigenvalue of 
of algebraic multiplicity equal to iV + 2. Therefore, since our structural assumptions 
ensure that srf v is a relatively compact perturbation of &f^> depending analytically on v 



(see Appendix B.l), there exist an analytic mapping v h-> I1(z/) where n(V) is a spectral 
projector for stf v of finite rank N + 2, and coincides with the orthogonal projector onto 
span($Q, . . . , <fr° N , <&° N+1 ) at v — 0. By Kato's perturbation method [HI pp. 99-100], 
we thus construct dual bases (<&q, <&i, . . . , & V N , & v N+ x) and *i> • ■ • > *jv> *at+i) °f> 
respectively, R(Il(z/)) and R(II(z/)*), which depend analytically on v in a real neighborhood 
of zero. The part of sf" on the finite dimensional subspace R(Il(z/)) is determined by the 
matrix 

£>":= ((^-^)) ( 
Similarly as .c/ 1 ', this matrix has an expansion 



0<a,fS<N+l 



D V = D^ + zifc^DW - k 2 u 2 D& + 



By using that 
^<f>° = 



cm, 



fc $o , 1<P<N, af {0) &% +1 = -c P k<f> 



o 

o • 



and (*°, $q) = 1, we get that 



= 0,1 < a < N, (^ (0) )**^ +1 = 



/0 



D (o) 



-/cc 



-kc MA 



-kcp i 
/ 



Using in addition Eq. (22), which equivalently reads 



(23) ^°)U fc + 



we see that for all a G {1, ... , iV}, 



*8> = o. 



hence 



V 



* 


* ... * 









* 








(Here above, we have also used that d u ^ u a ■ = for all a, and *° • ^W^,*" = 

for all a > 1.) Moreover, we claim that the upper-left entry of Z^ 1 ) is 
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Indeed, this equality comes from (23), and the only other nontrivial term in the upper- left 
entry of is 

i(*°-^)^($;) M >, 

Since (g/^)* (ty®) belongs to span(*°, . . . , ty° N , ty° N+1 ), that term will cancel out provided 
that for all a G {1, . . . , N + 1}, 

(24) (*° • d v (& ) w=0 ) = 0. 

This we can arrange, up to a harmless modification of $g. Let us explain how. Using 
that = 0, we see by expanding 

that 

n(o)K (1) $° + 4^ (0, ^(*S)i-=o) = ^ (1)$ o + ^ (0) ^(^) k=0 , 



or, using again ([23]) and that <3>$ G R(I1(0)), 

n(o)(-^(°)u, + I^(°)^(^)|, =0 ) = -^ (0) U fc + l^)^^)^. 

%K IK 

This shows that 

*f® (dv(*o)\u=o ~ ik\J k ) G R(n(0)) = ker((^°>) 2 ) = ker((^°)) 3 ) , 

hence 

0„(*5)i„=o - ifcu fc eR(n(o)). 

This means that there exist numbers z a such that 

d u {<f> u Q ) w=0 - ik\J k = z a $° a . 

If we substitute 



for $q, we thus have that 



^($ )| y= o = ik\J k 
(In order to keep duality we also substitute 



for ^> v a with a G {l,...,iV + 1}.) Forgetting the tilda, this implies (24) for all a G 
{1, . . . , N + 1} because 

«-U fe ) = 0. 
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As to the other diagonal block in it reduces to 

((*°-^ (1) ^))i<^<iv + i 

because (^)*9° a = for a > 1. It remains to compute the first column, starting from 
second row, in D^ 2 \ At a-th row we find 

1 «9,(<^V=o> + I (d u (* u a ) W =oi^ (1) K + ^ ^ (0) ^(*o)|,=o)) + <*° "^ (2) *§> 



where in fact the middle term is zero because of (24) and the fact that 

(&(*£W*S> = - «-^(*o)|,=o) 
by duality. Collecting together the results of the above computations, we get that 



D v :-- 



V 



tends to iK ) : = 

i) u fc + ^9 fl U)i) 



!) 9 e U)^) 



V ((j-^i-kwu. + ^^u)) 



-fcCMi ■ ■ ■ - kCM N 
(((^ (1) Um,) q )) l<a,/3<N 



((j-^.^Um;)!. 



—kcp 



(K (1) u P ) w ) 



((J-^-^'Up) / 



g/;i<P<N 

when z/ goes to zero. We are now going to check that the matrix of the modulated system 
Q linearized about (k, M = (U), P = (|U- J _1 U)) is nothing but 5 (0) - cIjv+2- Since 

1 S(z/)- 1 D i/ S(z/) 



with 



£(V) 



ikv 

( \ 





V 







/ 



the existence of a non-real eigenvalue of .D^ — cIat + 2 would imply the existence of an 
eigenvalue of D u bifurcating from zero into the right-half plane. Before linearizing it, let 



write (21) in the simplest abstract form 

drk = — dx(ck) , 
!2.V; { d T M = d x (JE^ k \V}), 

d T p = d x (u -Ejr( fe )[u] + ije^[v}} 
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where 

^T (fc) (U, <%U) := JT(U, £;<9 e U) . 



In quasilinear form, (25) reads 
(26) 

drk = —cdxk — kckdxk — kcugdxMp — kcpd x P, 
d T M = (3E^[V]) k d x k + (JEJf^\U\) Mp d x M p + (JE^[U]) P d x P , 
d T P = (U ■ E^( fc )[U] + L^ k \\J}) k d x k + (U • EJ^^[U] + LjeW[U]) Me d x M p 
+ (U-E^( fc )[U] + lJ^^{\J})pd x P . 

The right-hand side in the first row here above is indeed the first component of 

(D<® - cl N+2 )d x (k,M 1 ,...,M N ,P) T . 

It remains to identify the other rows. If a denotes any one of the variables M a or P, we 
have 

(JE^T (fc) [U]) a = (JHess^T (fc) [U]U a ). 



Recalling that 
and observing that 
we thus find that 



s/m = j(a(°) + kd e AP>) 

A (0) = Hess^T {fc) + cJ _1 , 
(JEJT«[U]) a = (^ (1) U a ) - c(U) a . 



Besides, we have 

(JE^( fc) [U]) fe = (JHess^ fc )[U]U fc ) + (JA«d U) = (^ (1) U fc ) + K (2) <%U) 
where we have also used that 

stfW = j( A W + fc&A^) . 



This shows that the second row in (25) equivalently reads in quasilinear form 

<9 T M = -cd x M + (^ {1) U fc + ^ 2) d e \J)d x k + (^ (1) U M(9 > d x M p + (^/ (1) U P ) d x P , 

in which the right-hand side clearly coincides with the a-th components, a £ {1, ... , iV}, 
of 

- cI A r +2 )9x(A;,M 1 ,... ) M JV) P) T . 
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In order to check the last component, we first compute that 
(U-EjrW[U] + UrW[U]) a = (U a ■ E^( fc )[U] + U ■ Hess^ fc )[U]U a ) 

+ (^^°>™«) "(!>") 

= (U ■ (J-V«U 4 )) - c (U • J- x U a ) - (U ■ kd e (A^V a )) 

+ <ra: (w ° )(w 4 

where again all derivatives of J$? are evaluated at (U, kdgXJ), and in fact the last four 
terms cancel out by definition of Therefore, using the symmetry of J, we obtain that 

(U-E^ fe )[U] + LJT«[U]) M , = ((rtj.^U^) 

and 

(U-E^T (fc) [U] + L^ (fc) [U]) P = - c+ ((J _1 U) -^ (1) Up). 

Concerning the derivative (U ■ EJ^f^[U] + LJ^^^U]}^, the computation is similar, with 
a few more terms: 

(U-EJTWfU] + L^W[U]) fc = ((J- X U) -^ (1) U fe ) + (U-AW^U) 

+ <ra; (w " )(W ) 

= ((J^U) ■ (^WU fc + ^d e V)) 
after another integration by parts and by definition of (plus the symmetry of J). □ 

2.3 Small amplitude wave trains 

A natural question is whether at least small amplitude waves are stable. Our main purpose 
here is to extend to our general, abstract framework, the necessary stability conditions that 
were exhibited by Whitham [211 PP- 489-491,512-513], basically for scalar problems. Small 
amplitude wave trains correspond to a near-linear situation. More precisely, solutions to 
([5]) involving only small oscillations (with bounded wavenumber) around a mean value M 
are expected to be well approximated by solutions to the linearized system 

(27) d t V = J? (HessJ?[M\U) , 
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with 

d 2 W d 2 W ( d 2 d 2 W 

(Hess^[M]U) a = 4r%r U e + Jr Z U t>* - D * Jr mt U ? + J Zr U ^ 
dU a dUp oU a dUp, x \oU atX dUp oU aiX dUp iX 

where all second order derivatives of J4? are evaluated at (M, 0). 

Let us concentrate for a while on periodic wave solutions to (JS]) of small amplitude, 
that is of the form XJ(x,t) = XJ_(kx + cut) with U being a 1-periodic function, (U) = M, 
and ||U — M||loo small. If in addition we assume that Jt? is an even function of (which 
is always the case when the contribution of to Jf? comes from a kinetic energy), the 
equations governing U 

( dW ( dW \ d£> \ 

(28) d e — (U, kd e V) - kd e — (U, kd 6 V) + c — (U) = , a E {1, . . . , N} , 



\dU a "-> °~ J °\dU a ^-' °- J J dU t 

with c = —u/k, are obviously symmetric under 9 t— > —9. In this case, we may assume 
without loss of generality that U is an even function of 9. Then, denoting by a = 
2|| (Ucos(27r-))|| an approximate amplitude of the wave, successive Lyapunov-Schmidt 
reduction arguments show that U can be expanded for small values of a as 

(29) U(fc,M,a;0) = M + aU^Jfe.Mjfl) + a 2 U 2 (£;,M; 9) + a 3 V 3 (k,M;9) + o(a 3 ) 

with (UJ = for all m = 1, 2, 3, (Uj • KU m ) = for m = 2, 3, (U 2 ■ KU 3 ) = whatever 
the NxN matrix K (in fact, since U is even, U^fc, M; 9) = 2 cos(27r0) (U^fc, M; •) cos(2tt-)), 
and JJ_ 2 (k, M; 9) = 2cos(47r6 l )(U 2 (A;,M; ■) cos(47r-))), and the frequency uj can also be ex- 
panded as 

(30) u{k, M, a) = u {k, M) + a 2 u 2 {k, M) + o{a 2 ), 

where Co>o(&,M) is the frequency of periodic wave solutions to ( p7| ), determined by the 
(linear) dispersion relation 

det{kA{2i7ck) - u {k,M)J- 1 ) = 0, 

(A(2i7rfc)) a fl := — + i (2nk) — + (2tt£;) 2 , 

all derivatives of Jif being evaluated at (M, 0). Furthermore, u^ffc, M) can be expressed 
in terms of mean values involving U 1 and U 2 as 

, A/n . Wh ■ d e (5^[Mmi,U 2 ) + ^^[MKU^U^U,))) 

(31) u 2 (k,M):=-k ^ 3 _ ld&) , 

where 5 3 J^f^ and 5 4 J4? ( - k ^ denote respectively the third and fourth order variational deriva- 
tives of J4?( k ) : U 1 — y J4?(XJ, kdeXJ) (up to this point, we have preferred the notation Ejtf? 
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for the first derivative and HessJrf? for the second order one 8 2 J^). The actual deriva- 
tion of (31) follows the lines of computations made by Whitham in j2TJ pp. 472-475^ 
first for KdV and then for water waves in arbitrary depth, except that here we have no 



explicit formula for ujq. The reader may check that (31 ) is consistent with values found by 
Whitham, namely 2nu 2 (k,M) = -3c 2 /[32 (2nk)} for (gKdV) with f(v) = c a (\v 2 + \v 3 ) 
(see [211 PP- 463,473]), and also 2ttu; 2 (A;,M) = -3a (2nk) + 24a 2 M 2 /(2nk) for (gKdV) 
with f(y) = \v 2 + av A , the zero-mean periodic solutions of which are basically governed 
by the same equation as periodic solutions to the Klein-Gordon equation considered in 
PCE1 pp. 486-487]. 



Derivation of (31 ). By plugging (29 ) and (30 ) in (28 ), we get successively, by increasing 



order in powers of a, 
(kA(kd e ) - w (k, M) J- 1 ) d g TJi 
(kA(kd e ) - u (k, M) J- 1 ) <%U 2 
(kA(kd e ) - oj (k, M) J" 1 ) d e V 3 



dedkS^^iMm^V,)), 

d e (kS 3 J^^[M)(U 1 ,\l 2 ) + lk5 A J^^[M){V l ,V 1 ,'U 1 ))) 
+ co 2 (kM)^ 1 deU 1 . 



By taking the inner product of the last equation with de]J_ 1 , we see that the left hand 
side vanishes because of the first equation and of the self-adjointness of the operator 



(kA(kd e ) - u (k, M)J" 1 ), hence pi Note that for the very same reason, one can draw 
information on uo^{k, M ) from the first equation. Indeed, by differentiating it with respect 
to k and taking the inner product with de]J_ 1 , the term involving dg\J_ 1 k vanishes, and we 
find that 



(32) 



u 0;k {k,M) 



(deHi ' {kA(kde))kdeHi) 



More explicitly, the numerator here above reduces to 

(V^dOj^U!,^!)) - 3k 2 (V 2 Jx J^(M,O)(d e U 1 ,d 3 U 1 )) . 

Returning to a more general small amplitude modulated wave train, the previous 
remarks show that we can expand its lower order term U°, which we merely denote by U 
in what follows, as 



U(T,X, 



(U)(T,X) + a(T,X) Ui(fc(T,X),(U)(T,X);0) 

+ a(T, Xf U 2 (k(T, X), (U)(T, X); 9) 

+ a(T,X) 3 U 3 (k(T,X),(U)(T,X);9) + o(a(T,X) 3 ) 



Our aim is to show that, for small a, the modulated system (21) associated with U 
decouples into the lower order system 



(33) 



<9 T (U> = Jd x (VuJr«U},0)) 



2 We warn the reader that we have taken the opposite sign for ui compared to the one chosen by 
Whitham and normalized periods of profiles to one instead of 2ir. 
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and the 2x2 system (pointed out in 
(34) 

d T k -u 0>k (k, (U)) d x k 
d T (a 2 ) - Lo M k > ( u ))« 2 d xk 



p. 490] when (U) = M) 



-u 2 (k, (U» d x (a 2 ) 
-Uo, k (k, (U))d x (a 2 ) 



w ,m(A;,(U)) d x (U) 
0(a 2 )d x (U). 



As a consequence, a necessary condition for (plj) to be hyperbolic in the small ampli 



tude limit is that both (33) and the lef-hand side of (34) with frozen mean be hyper- 
bolic. Regarding the latter, this requires that u 2 (k, (XJ))u>o,kk(k, (U)) be positive, in 



which case the characteristic velocities of (34) are 0(a) perturbations (—oJo,k(k, (U)) ± 
a^ui 2 (k, (XJ))uJo,kk(k, (U))) of — uio,k(k, (U)). As to the former system (33), it is for in- 
stance nothing but the first order conservation law 



(35) 



d t {v) + d x p((v)) 







if the original equation is (gKdV), that is Eq. Of course (35) is always hyperbolic, as 
any (real) first order conservation law. However, the hyperbolicity of the reduced system 



(33) is not automatic in general. In particular, when we start from the Euler-Korteweg 
system, we receive as (33) the Euler system. Therefore, a necessary condition for the 



modulated system to be hyperbolic in the small amplitude limit is that the mean value 
of the wave train be a stable state of the Euler system. According to Theorem [TJ this 
implies that small amplitude periodic solutions whose mean value is an unstable state of 
the Euler system are themselves (spectrally) unstable. 



Derivation of (33) and (34). There is not much to do about the first equation in 



(34). Indeed, the first equation dxk + d x (ck) 
-u {k, (U» - a 2 u 2 (k, (U)) + o(a 2 ) yield 



in (21) together with the fact that 



ck 



d T k - (w ,fc + a 2 u 2 ,k)d x k - u 2 d x (a 2 ) = (w ,m + a 2 u 2M )d x (U) + o(a 2 



33 
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and the terms a 2 u) 2 ,k and a 2 cj2,m are negligible compared to the 0(a) perturbation of 
— wo,fc we are expecting (this was already pointed out in [2H p. 490]). So the main points 
consist in showing that the middle equations <9t(U) = J d x (G) in plj ) do reduce to 
when a goes to zero, and that together with the last equation 8t{Q) = d x (S) in 
they simplify into an equation for a 2 whose principal part amounts to 

d T (a 2 ) - d x (u 4k,(V))a 2 ) = 0(a 2 )d x (V). 

Let us start by expanding G in powers of a. Note that we shall need an expansion not 
only for (G) but also for (U ■ G), which arises in S. Recall indeed that 



G = 5J^ {k) [U] , S = U • G - JT(U, kd e U) + {kd e U a 



dU nr . 



(U, WJ). 



We thus see that 

G = dje^iiV^+ad^^liV^V^ + a 2 (<^ (fe) [(U)](U 2 ) + itflr (fc) [(U)](Ui, ILJ) + o( 
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with 

5^[(v}} = Vujr((u),o), 

and 

(<far«[(U)](Ui)) = 0, (5 2 je^[(V)](V 1 )) = 0, 

since 5 2 J^ fe )[(U)] is a constant-coefficient differential operator, and Ui, U2, as well of 
course as their derivatives, have zero mean values. Therefore, it just remains 

(G) = Vu^((U),0) + |a 2 (^ fc )[<U)](Ui,Ui)) + o(a 2 ) . 



By neglecting all the 0(a ) in <9r(U) = 3 dx(G), we thus receive (33). Let us turn to 
the expansions of (Q) and (S). The expansion of (Q) is readily seen to be given by 

(Q) = i(U) ■ 3- 1 (U) + i a 2 (Ua • J" 1 ^) + o(a 2 ) , 

and we also find that 

(U ■ G) = (U) ■ (G) + a 2 (Ui • ^^[(UJKUO) + o(a 2 ) . 

(We have used again that Ui and U2 have zero mean values.) Furthermore, we have 

(JT«[U]) = Jt?M[(V)] + §a 2 (<^*>KU)](Ui,Ui)) + o(a 2 ), 

and finally 

f) yf> f) %P ft ffi> 

(kd e U a )—(U,kd e XJ) = a(kd e U la )—((XJ),0) + a 2 (kd e U 2a ) — «U>, 0) 

+q2 (^ (^S ( ^ o)+w & ((d,,o) ) + o(a2) ' 

hence 

(kdeU a )-j^(\J,kd e V)\ = a 2 k 2 (Vl^mMdeUiideUi)) + o(a 2 ) . 

Now, taking the inner product with (U) of J -1 <9t(U) = dx{G) and subtracting to 
9t(Q) — dx(S) we see that 

d T (| a 2 <U X ■ J- 1 ^}) = 9x((U-G) - (^ {fe) [U]) + (^kdoUJ^tUtkdoV) 

-(U)-dx(G) + o(a 2 ). 

As expected, the zeroth order terms in the right-hand side cancel out, and we receive after 
several simplifying operations 

d T (a 2 (Ui • J-'Ui)) = 0x(a 2 (V^((U),O)(Ui,Ui) 

+ 3 k 2 V 2 J ^((U), 0)(a e Ux, Wi)) 
+ a 2 (5 3 ^ fc )[(U)](U 1 ,U 1 )) 9 X <U> + o(a 2 ). 
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In order to conclude, let us recall (32) hence, since the dependency of U 1 is a cosine 
function, 

(Vl^((\J)MViiVi) + 3fc 2 V^((U),0)p 1 ,() 9 U 1 )) 
u 0jk (k,M) = — — — in jt r • 

Therefore, the equation above reads 

d T (a 2 (V 1 -J- 1 V 1 )) + d x (co 0ik (k,M)a 2 (U 1 -3- 1 V 1 )) = 0(a 2 )d x (U) + o(a 2 ) . 
The factor (Ui ■ J _1 Ui) can now be eliminated by using the equation on k. 

3 Application to the Euler— Korteweg system 

In this section, we concentrate on the Euler-Korteweg equations, (JTJ in Eulerian co- 
ordinates, and ([2]) in mass Lagrangian coordinates. We derive Whitham's modulation 
equations for both systems, and show that away from vacuum, the modulated systems 
are equivalent through a mass Lagrangian change of coordinates. Thus, it is sufficient to 
check the hyperbolicity of modulation equations for either one of these systems in order 
to determine whether our necessary condition for the spectral stability of periodic waves 
under small wave number perturbations is satisfied. 

3.1 Periodic traveling waves 

Periodic traveling wave solutions to (jlj and (j2]) are respectively of the form (p, u) = 
(R, U)(x — at) and (v, u) = (V, W)(y + jt), with a one-to-one correspondance between the 
two frameworks encoded by 

R(OV(Z(0) = 1, U(0 = W(Z(0)^ = R = Y(zy 

Up to translations, these periodic traveling waves generically arise as four-parameter fam- 
ilies. Natural parameters are 

• their speed, that is a in Eulerian coordinates, and —j in Lagrangian coordinates, 

• a first constant of integration, which turns out to be j in Eulerian coordinates, and 
u in Lagrangian coordinates. 

• two other constants of integration/Lagrange multipliers, which we denote by A and 
li, in the profile equations. 



To be more precise, the profile equations read 
(36) 

L V-p-* 1 - ) \ 1Xi i ±L £.) — ~ " 
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R(U - a) = j , ( W - jV = a, 

(L p J?)(R,Rs) = - A i {LJ)(V,V C ) = - fx, 



where the Lagrangians Jzf and I are defined by 



(37) # := £-£- liP) t:=£ -^.-Xv, 

and 

L % ■= o — — - & L £ - v — - £ 

are (obvious) first integrals of the Euler-Lagrange equations E p J?f = and E v £ = 
respectively. In addition, there is a simple relationship between the mean values of Eule- 
rian profiles and of Lagrangian profiles. Indeed, if S is the period of a traveling wave in 
Eulerian coordinates, the period of its counterpart in Lagrangian coordinates is Z(E) (if 
Z is chosen so that Z(0) = 0), and we have 

1 r Z(E) 1 f z{B) - 

(R) := ^ / R(0 d£ = , (V):=— V(C) dC 



^(s)y ^ ' z(s)' 

(U):=±£u(0dZ = (^):=^y^V(C)dC = JL(i2tf>, 

hence the remarkable identities 

( R )-Jyy ^-^yT' W -7i?T- 

Note that these relations are of course compatible with the profile equations R(U — a) = 
j , W — j V = er. We warn the reader that for convenience we denote by the same brackets 
(•) the mean values with respect to £ and those with respect to (. This should not be 
too much confusing since we use different notations for the Eulerian and the Lagrangian 
dependent variables. 

3.2 Whitham's modulated equations 

As done previously in the abstract case, we look for solutions of Q and ^ having 
asymptotic expansions of the form 

(p,u)(t,x) = (po,uo)(st,ex,^(et,ex)/e) + e(px,ui)(et,ex,^(et,ex)/e,e) + o(e) , 

(v,u)(s,y) = (v ,w )(es,£y,(p(es,ey)/e) + e (vi,wi)(es,ey,<f>(es,ey)/e,e) + o(e) , 

where the profiles (po,Uo), (pi,u\), (vo,wq), and (vi,wi) are 1-periodic in their third 
variable 9. Denoting by T, X, S, and Y the rescaled variables (et, ex, es, and sy), we 
introduce the further notations 

n 

K : = $ x , := $t, a- := - — , 

A 
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k := <p Y , u ■= 4>s, j ■= -r- 
The existence of smooth phases $ and <fi requires the standard compatibility conditions 

(38) d T K + d x {(iK) = 0, d s k - dy(jk) = . 

We obtain equations for the leading profiles in both Eulerian and Lagrangian coordinates 
by plugging the asymptotic expansions in Q and Q, using that 

d t = ed T + ^ d 9 , d x = ed x + K d 9 , d s = ed s + ui d e , d y = ed Y + k d e , 

and retaining only the leading order terms. As expected, we recover the traveling profile 



equations (36)-(37), up to the rescaling that makes their period equal to one. In Eulerian 



coordinates, this gives 

(39) d e (p u ) -crdepo = 0, (u - a) d e u + d e g = 

where 



dS ( dS \ 

9o ■= -7^(po,Kd e p ) - KD e (—(po^depo)) , 

while in Lagrangian coordinates, 

(40) d e wo - j d e v = , j d e w + d e p = , 

where 

Po ■= ~ q^{vo, kd e v ) + kD e \^—(v ,kdgVo)J . 

The next order leads to Whitham's modulated equations. Indeed, collecting all the 
terms of order one in (fij) and (|2]) we get 
(41) 

-oKdepi + Kdeipoux + p 1 u ) + d T po + d x (p u ) = 0, 
-aKdeUi + Kd e {uQUij + K d e (A pi + B ) + d T u + u d x u + d x go = 0, 



(42) 
where 



j k d e v l - kd e w 1 + d s v - d Y w = , 
jkdgW! - kd e (a v 1 + b ) + d s u + d Y p = 0. 



d 2 S 

A := -KDe^oKDe + T , JT := ^(po, Kd e p ) , 
To := -^r{po,Kd e p ) - KD g ( (p ,Kd e p ) , 



— {Po ,KcJePo)-KD e 

d pdp x — mpo)dxp °- dx {w, 



d 2 <§ ( d<S \ 

B o := a a (Po,Kd 9 p )d x p - d x I -—(p ,Kd 9 po) - K D e (Jtf d x p ) , 
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d 2 e 

a := -kD e K kD e + 70, k := — (v ,kd e v ) 



d 2 e ( d 2 & 

~» : = -^( v o,kd e v ) - kD e I Q Q (vo,kd e v ) 



dv" z 



dvdv,, 



d 2 , 



dvdv, 



(v 0j kdev )d Y vo - d Y ( S^-(v , kdev Q ) ] - IcD 

\dv y J 



e {k OyVq) 



Integrating (41) and (42) with respect to 6 over [0, 1], we get the modulated equations 

d T (po) + dx(pouo) = 0, 



(43) 



d T (u ) + d x {\u 2 ) + d x (g ) = 0, 



(44) 



d s (v ) - d Y (w ) = 0. 
ds{w ) + d Y (p ) = 0. 



(Observe that the terms B and bo, the only ones involving cross derivatives of the leading 
profiles, play no role at this level.) Now, taking into account the compatibility condi- 



tions in (38), we only need to find a fourth averaged equation to have a complete set of 



modulated equations. We can proceed in two ways. 

As mentioned in the abstract case, the fastest way consists in averaging additional 
conservation laws satisfied (at least formally) by solutions of Q and Two possible 
choices are the conservation law of total energy (associated with invariance with respect 
to time translations, via Noether's theorem), and the (local) conservation law of Ben- 
jamin's impulse (associated with invariance with respect to spatial translations). These 
conservation laws read 



(45) 



d t {\pu 2 + S) + d x (^pu 3 + puE p S + d x {pu) 



dp x 







(46) d t (pu) + d x {pu 2 + P E P S + L P 

in the Eulerian framework, and 



0. 



(47) 



d s (^u + e) - d y [ uE v £ + (d y u) 



dv,, 



(4£ 



d s {yu) — d y (\u 2 + v E v £ + L v £) = . 



in the Lagrangian framework. Upon plugging the asymptotic expansions and averaging 
we get 

(49) d T {\p Q u 2 + Sq) + d x {\poul + p u g + K d e (p u ) —(p , Kd e p )j = 0, 
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(50) d T (pou ) + d x (poul + p g + K (d e p ) — (p ,Kd e p ) - S{p 0l Kd e po) 



dp x 



in the Eulerian framework, and 



(51) 



9s(l wl + £q) + dy ( w po - k (dgw ) q^-( v 0' kde v o) ) = , 



(52) d s {v w ) + d 



Y 



v po + £{vo,kdgv ) - k(dgv ) —(v ,kdgv ) J 







in the Lagrangian one. For simplicity, we have denoted by Sq and e§ the energies evaluated 
at (p ? Kdgpo) and (vo,kdgvo) respectively. At first glance it may look like we have five 
modulated equations in each framework, namely the compatibility condition in (38), the 
two equations in either (43) or (44), as well as (49)-(50) or (51)-(52). In fact, the profile 
equations in (39) imply that the averaged energy equation (49) is a consequence of (38) 
and Q g, and similarly Q-(|44|)-([52]) imply |>T]). 

A more indirect way to derive a fourth modulated equation consists in taking the inner 
product of (41) and (42) with (u ,po) and (u>o,t>o) respectively. Accordingly with the 
abstract case considered in Section [2j those choices are dictated by Benjamin's impulses 
pu and vu, of which the variational derivatives are respectively (u, p)* and (u, v) f , see 
Appendix [A] for more details. 

To summarize, we have the following. 

Proposition 1. Whitham's modulated equations associated with ([T]) and ^ read, respec- 
tively, 

• in the Eulerian framework 

' d T K + dx(o-K) = 0, 
9t(po) + dx(pouo) = , 
d T (u ) + d x (\ul) + d x (go) = 0, 



(53) <^ 



dS 

d T (pou Q ) + d x ( poul + p g + K (d e p Q ) —(p , Kd e p ) - S Q 



0. 



which is endowed with the additional conservation law 
(54) 



/ dS 
d T {\pvul + £ Q ) + d x (\pQul + p u go + K d e (p u ) —(p , Kd e p )) = 

and in the Lagrangian framework 

( d S k - dy(jk) = 0, 

<9s(^o) - <9y(w ) = 0, 
( 55 ) { d s (w ) + d Y (p ) = 0, 



ds(v w ) + d Y (- \ wl + vopo + £ 



de \ 
k (dgv ) g^-( v o, kdgVo)) 
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which is endowed with the additional conservation law 



(56) 



de 



dsihwl + Eo) + dylwoPo - k{d e w Q ) —{vo.kdevo) 



0. 



Now, we may go further and make the link between the Eulerian (53) and the La- 
grangian (55) modulated equations. Interestingly, even though it is all but a surprise, 
(55) can be viewed as a Lagrangian reformulation of (53). More precisely, we are going 
to show the following. 

Theorem 2. Away from vacuum, there is a mass Lagrangian system of coordinates chang- 
ing System (53) into (55). In particular, these systems are simultaneously hyperbolic. A 
sufficient condition for the hyperbolicity of (53) and (55) is that 

e := (£o) + | (wl) - | (w ) 2 
be a strictly convex function of (v, k, A/k), or equivalently that 

E:=(p )e = (4) 



1 / 2\ 1 (P0U ) 2 

2^ - 2~wr 



be a strictly convex function of (g,K,D/K), where 

v := (v ) , A := (v w ) - (v ) (w ) , 

Q-=(Po), D ■= (po)(u ) - (pqU ) . 

Proof. Let us recall that (R, U) is a l/i^-periodic solution of ( |64| ) if and only if (V, W) is 
a 1 //c-periodic solution of (65), with 



R(C)V{Z{0) = 1, U({) = W{Z(Z)), 



dZ 



R 



1 

k 



This implies in particular that 

K 



(57) 



(Po) 



k 



k 
K 



(vo) 



(Po) 



( w o) 



v{zy 

(pouo) 



K 



(Po) 



These observations enable us to make the relationship between (53) and (55) in the same 
way as it is usually done between the fluid equations in Eulerian coordinates and those 
in mass Lagrangian coordinates. As a matter of fact, the second equation in (53) states 
that (p )dX — (p uo)dT is a closed differential form, and thus an exact form in a simply 
connected domain. We thus introduce the 'rescaled mass Lagrangian coordinate' Y defined 
(up to a constant) by 

dY = (po)dX - ( Po u )dT. 
Setting S = T, this equivalently reads thanks to the last two relations in (57), 

dX = (v )dY + (w )dS, 
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hence the second equation in (55). We can proceed similarly with the other equations. 
The first one in (53) states that 



K dX — oK dT = (v )KdY + K((w )-a)dS 



is an exact differential form. Using the second relation in (|57)) and that (wq) — a = j (vq) 



we get that kdY + jkdS is exact, hence the first equation in (55). As regards the other 



equations, the third one in (p3l gives the fourth one in (55 ), and the fourth one in (53 ) gives 



the third one in (55) (this interplay between momentum and velocity equations is well- 
known for standard fluids motion, in which the conservation law of the momentum pu in 
Eulerian coordinates is associated with a conservation law for the velocity in Lagrangian 
coordinates, the other way round going from a conservation law for vu in Lagrangian 
coordinates to a conservation law for u in Eulerian coordinates being less classical but 
still true, as far as smooth solutions are concerned). In order to justify the correspondence, 
it is convenient to rewrite these equations in a simpler way. This is done by using the 
profile equations, which give 



Po(«o -cr) = j , g 



j 2 OS 
V> — 77-2 > K (d e p ) —{p ,Kd e p ) - & 
2po dpx 



0, 



2p 



ppo - A 



w - jv = cr , po = -A - j D , Bo - k (d e v ) i^—(v , kd e v ) 

8Vy 



Ifvl + Xv + p, 



so that the last two equations in (53) and (55) respectively read 



(5* 



d T (u ) + d x (p - \ o 2 + a(u )) = , 

' 2 > — A — j 2 ( — 



d T (pou ) + d x (poUq) - A - y 



(59) 



d s (w ) - d Y (X + j> » = 0, 

ds(v wo) + d Y {p - \a 2 - j a (v ) - j 2 (ug)) = 0. 



The first equation in (58) means that 



(u )dX + (la 2 - p - a(u )) dT = (v w )dY + (ja 2 - p + ((w ) - a) dS 
= (v w )dY + (\a 2 - p + ja(v ) + f (v 2 )) dS 



is exact, which is equivalent to the second conservation law in (59). The second equation 



in (58) means that 



( Po uo)dX + (- (p u 2 ) + A + j 2 (J~yj dT 
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(w ) dY + 



- (p u 2 ) + A + j 2 (iXj dS=(w )dY + (A + j 2 (v ) 2 ) dS 



is exact, which is equivalent to the first conservation law in (p8|). This finishes to prove 



the equivalence between the two modulated systems (53) and (55) (as long as (p ) and 
(vq) do not vanish). As a consequence, these two first order systems of conservation laws 
are simultaneously hyperbolic. A sufficient condition for the hyperbolicity of (55) has 
been pointed out by Gavrilyuk and Serre p2] in terms of the average energy 



e := (£q) + \ (w\) - \ (w ) 2 
which satisfies the generalized Gibbs relation 
(60) de = - pdv + 6 dA; + j dA 

with p := (po), v := (v ), and 

de 

[V , k(j e V 0/ 



(61) 



6 := ( {d e v ) q^-( v o, kdgv ) ) , A := (v w ) - (v ) (w ) . 



Namely, it is shown in [12j (alternatively, the reader might take a look at [3]) that, if e 
is a strictly convex function of (v,k,A/k) then (55) is Godunov-symmetrizable and thus 
hyperbolic. We may note in addition that, by a standard argument on convex functions 
(viewed as supremum envelopes of affine functions), the convexity of e as a function of 
(v, k, A/k) is equivalent to the convexity of (p ) e as a function of ((po), k(po), (po)A/k). 
An easy calculation (using relations in (57) and similar ones) shows that 

1 (Po^o) 2 



(p > e = (Sq) + -(p «o) - 2 / v 



E, k(po) = K, (po) 



A 



K 



((po)(u )-(poU )) 



This shows that a sufficient condition for the hyperbolicity of (53) is the strict convexity 
of E as a function of (g, K, D/K) where g := (po), D := (po)(uo) — (po^o)- Of course our 
generalised Gibbs relation has its counterpart in terms of E. It reads 



(62) 

where 

(63) 



dE = gd£ + 9 dK + - dD 

Q 



(go) + 



r 



i 

? 



i 

7 



Q 2 



Equation (62) can be derived from (60) as follows. Using that 

E = £e, v = 1/g , K = kg , D = g 2 A . 



we readily get that 
dE 



e + pv - 



KQ 



jD 



dg + 6 dK + - dD . 

Q 
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It just remains to express the factor of dg in terms of 'Eulerian' mean values. Indeed, 
we can see the 'temperature' either as a 'Lagrangian' mean value (by definition, see 
Eq. (61)) or as an 'Eulerian' one: it turns out that 

6 = ((d 9 po)—(p ,Kd e po)) . 

As a matter fact, by the Eulerian profile equation we have 

I dS \ I j 2 \ I j 2 v 2 
K UpePo) Q^-{po,Kdepo)\ = ( S - — - pp - A j = g ( e ^ - p - Xv 

by the relationship between Eulerian mean values and Lagrangian mean values we already 
used several times, hence by the Lagrangian profile equation 

K ((depo) ^-(po,Kdepo)^ = gk ^(d e v ) -^-(v , kd e v )^ = KB. 

We can now work out the factor of dg by observing that 

jA = (w j(v -v)} = (w (w - (two))) = (w%) ~ (w ) 2 , 
so that by definition 

e = ( £ o) + -y • 

In addition we have that 

p v = (p v ) + j A 

since 

(PoVo) = ((p-j 2 (^o-v))w ) = pv-(j 2 K-v) 2 ) = pv-((w -(«)o» 2 ) = Pv-j'A. 
Therefore, 

e + p v = {£q + p d ) + -y ' 

and by our usual trick, 

(£Q + po v ) = - (So + po) , 
Q 

where po has to be expressed in Eulerian coordinates in the right-hand side. This amounts 
to writing 

(<?o+Po) = U, - X - 3 -\ = KQ- 3 ^- (-\+pg= KQ-^ (-)+Q (go + A 
\ Po/ 2 \po/ 2 \p Q / \ 2p z 



using once more the profile equations. So we have 
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and noting that 



Q 



1 

Po 



(fv 2 ) = j 2 v 2 + j A 



we eventually obtain the claimed formula (62) with g defined by (63). 



□ 



Still, it is not obvious at this stage that (53) or (55) are really evolution systems in 
closed form. We shall observe on numerical experiments that evolutionarity may indeed 
fail. For the moment, let us just note that if evolutionarity happens to fail, it does so 



'simultaneously' for (53) and (J55J). As a matter of fact, at fixed (T, X), we know that 

(p ,u )(T,X,9) = (R,U)(T,X,e/K(T,X)) 
where the profile (R, U) is a 1/X-periodic solution to 
R(U -a)= j, 



(64) 



Rz—(R,Rt) - <?(R, R$) + ^— + fiR = - X 
OPx 



Here above, the sign = means equalities for functions of (T, X) only. The quadruplet 
(j,o~,\,fi) is a natural set of parameters in (64), which we expect to determine properly 
the wave number K as well as the l/f^-periodic solution (R, U) to (64), up to translations, 
hence also all mean values involved in (53) (54). Similarly, 

(y , w )(S, Y, 6) = (V, W)(S, Y, 9/k(S, Y)) 



and all mean values involved in (55 ) ( |56| ) are expected to be fully determined by (j, er, A, (j) 
through the profile equations 



(65) 




(V, V ( ) - *(V, V C ) + 



J 2 V 2 



XV 



LI. 



By the one-to-one correspondence we have pointed out between Eulerian and Lagrangian 
periodic orbits, the mapping (j,a,\,Li) H- (k, {V), (W), (VW)) will be locally invertible 
if and only if (j, a, A,/i) i-> (K, (R), (U), (RU)) is so. Again, this is not exactly always the 
case, as we shall see on specific examples in the next section. 



4 Nature of modulated equations in practice 

Our purpose here is to investigate the hyperbolicity of Whitham's modulated equations 
for the Euler-Korteweg system. First of all, let us point out that the parameter a does 
not play any role in that matter. This is due to Galilean invariance of the Euler-Korteweg 
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system, unsurprisingly. To be more precise, let us rewrite the modulated equations in nat- 
ural coordinates from the 'thermodynamical' point of view. In the Lagrangian framework 
they read 



(66) 



d s k - dy(jk) = 0, 
d$v — dy w = , 
<9 5 w + <9 y p = , 
<9 5 (A + vw) + dy (e + vp 



kQ 



±w 2 

2 v 



together with the generalized Gibbs relation (60) de 



2jA) = 0, 
pdv + 9 dA; + j dA. We easily 



check that (66) is invariant by translations of the form 



(S, Y, k, v, w, A) !->■ (S, Y, k, v, w - a, A) 



for any (constant) velocity q_. Substituting o + jv for w, we thus see that the equations 
in (66) are unchanged if a is replaced by a — q_. Therefore, the eigenvalues of the lin- 
earized equations about a reference 'state' (k, y, er, A), and those about the translated 
one (k, y, 0, A) coincide. A similar argument works in the Eulerian framework too, see 
Appendix [C] for more details. 

From now on, we concentrate on the Euler-Korteweg system in mass Lagrangian 
coordinates ([2]), with an energy of the form (|4]). We recall from (65) that a periodic 
traveling wave solution (v,u) = (V, W)(y + jt) to ([2])-(|4]) is characterized by a periodic 
profile (V, W) that must be solution of 

(67) W — jV — a, ^(V)V ( 2 + \fV 2 + XV - f(V) = -p, 

where a, A, and /i are constant of integrations. 



Remark 1. The phase portrait of the ODE on V in (67) does not depend on the specific 



capillarity function v i— >■ k(v ), provided that it stays positive. In fact, up to a rescaling in 
(, that ODE reduces to 



iy 2 

2 ' 



f(V) - W 



XV 



Remarkably enough, this equation is also the integrated profile equation for traveling wave 
solutions of speed —j 2 to the generalized Korteweg-de Vries equation (gKdV) 

d t v + d x p(v) = —d^v . 

(This relationship between the traveling waves of the Euler-Korteweg equations in La- 
grangian coordinates and those of the generalized Korteweg-de Vries equation has been 
known for a long time, and is investigated in more details for instance in JT3]/ J However, 
when one turns to specific examples for p in the generalized Korteweg-de Vries equation, 
it is most often to consider power laws p = t> 7 . By contrast, we consider here laws that 
go to infinity at zero, or more generally at some co-volume b, and to zero at infinity. 
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Using in particular the first equation in ((67) above to substitute a + jV for W in (|55|) 



and reformulating the last two equations in (55) as in (59), we can write the modulated 
equations as 

(68) d s k-d Y (jk) = 0, 

(69) d 3 (V)-dy(j{V)+(r) = 0, 

(70) d s (j(V) + a) - d Y (\ + f(V)) = 0, 

(71) d s {a(V) + j(V 2 )) - d Y (W + ja(V) + f(V 2 ) - /i) = 0. 



This is just an alternative formulation of (66) in terms of the 'natural' parameters 
(cr,j, A, /i). Our main purpose is to investigate the hyperbolicity of (68 71 ), which by The- 
orem [I] is a necessary condition for stability of the periodic wave, provided that nearby 
periodic waves be parametrized by {k, (V),j(V) + a, u(V) + j(V 2 )). We can hardly check 
these properties - evolutionarity and hyperbolicity - analytically, since neither the wave 
number k nor the wave profile V is known explicitly in terms of (a, j, A,/i). However, it 
is not difficult to check them numerically. To make numerical computations more trans- 
parent, we are going to change (A,/i) for more convenient parameters, under suitable 
assumptions on the pressure law v t— > p(v) = —f'(v). In this respect, we shall start with 
shallow water type pressure laws, and turn to the more complicated, Van der Waals type 
pressure laws afterward. We shall see that the form of the capillarity n(v) also plays a 
role - if a rescaling along orbits does not change the phase portrait, it does change mean 
values. Our minimal assumption will be that k : (b, +oo) — > (0, +oo) is a smooth function 
for some nonnegative 'co- volume' b. 

Before going into specific examples, let us say a few words about two asymptotic limits, 
namely the small amplitude regime, and the solitary wave limit. The former has been 



analyzed in some detail in |2.3[ In particular, it has been pointed out that a necessary 
condition for the hyperbolicity of modulated equations about a small amplitude periodic 
wave for the Euler-Korteweg system is that the Euler equations be hyperbolic at the 
mean value of this wave. This is why we shall not try to get too close to center points 
of the wave profile equations in the numerical computations that follow: we readily know 
that the small amplitude wave trains about those points where the Euler equations are 
not hyperbolic are unstable. (As to the other necessary condition, namely the positivity 
of W2Wo,fcfc, we leave it aside for the moment.) The solitary wave limit corresponds to 
when the wavenumber k tends to zero. As noticed in earlier work (see for instance [8j [9]) 
for which explicit computations can be made - involving elliptic integrals - , we expect 
modulational instability for waves of small wavenumber when the endpoint of the limiting 
solitary wave is an unstable state of the Euler equations. Thus we shall not try, in 
numerical computations, to get too close to solitary wave orbits either. 
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4.1 Shallow water type pressure laws 



4.1.1 Parametrization of periodic waves 

Here we assume that p : (5, +oo) — » (0, +oo) is smooth and strictly convex, with 

limp(t>) = +oo , lim p{v) = , 

hence in particular p is monotonically decaying to zero at infinity. The shallow water case 
p(v) = 1/v 2 is the main application we have in mind. If is to denote the endpoint of 
a solitary wave with j ^ and A as constants of integration, it must satisfy 

f < -p'iVoo) , A = -j 2 Voo - p(Voo) , 

in which case there is exactly one other solution t> G (foo, +oo) to A = —j 2 v — p{v) . Now, 
inside the homoclinic loop connecting Voo to itself, there is a collection of periodic orbits 
surrounding vo, which are determined for instance by their trough v* (see Figure 4.1.1). 
More precisely, if G (foo^o) there is a unique periodic orbit passing through v* and 
solving the ODE 

\k{V)V 2 + IfV 2 + XV- f(V) = -fi, /x := f(v m ) - \M - \v, . 



P=-f(v) 



-X 





/ slope - 


-P 




-M 




v' 













Figure 1: Phase portrait for convex pressure law 
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Let us consider the mapping 

$:OcR 3 ^R 3 
(j>oo,v*) ^ (j, A = -fv^ -p(«oo), At = /(«*) - |j 2 ^ - Au*), 

with 

^ = {0') ^oo, «*) | < j 2 < -p'(foo), 6 < ^oo < < ^0 ; fVoo + P(Voo) = fv + p{v ) } ■ 

It is a diffeomorphism onto the open subset A of IR 3 made of parameters (j,X,fi) for 
which we do have a periodic wave. From now on, we parametrize periodic waves by 
(ji v oa,v*) G fi instead of (j, A,//) G A. Given (j, u*) G and (j, A,/i) = $(j, -Woo, t>*) 
there is a unique t> * (the peak of the periodic wave) greater than t>* such that 

/K) - §iV) 2 -\v* = ^ = /(«.) - ±a. 2 - a^ . 

The wave number of the corresponding periodic wave is given by 



i _ r l «<«> iv 



2k I, V 2(/(«) - ijV - At, - p) 
and mean values can be computed by quadrature in a similar way. In particular, we have 



IT 



2(f(v) - IjhP-Xv-fi) 



^ = 2 *> ^/ 2(/ (,)-^-A,-,) d - 

An elegant way to remove (integrable) singularities at endpoints in the integrals above 
has been pointed out in |5J. Indeed, factorizing the denominator as 

f(v) - W^ 1 -Xv-fj, = (v* -v)(v- v„) ip{v; j, Woo, «*) , 

and introducing the change of variables 

~/s m ■ \ v* + v* v* -v* 

T : (0; j, Voo, u*J i ^ f •— — 1 — sin 9 , 

we find that 



i r' 2 j K(-r{8;j,voo,v*)) d0 



2k J-n/2 V 2^(f(6 l ;j> 00 ,^);j,w 00 ,^) 
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(V 2 ) = 2k 



it/2 



n(i / (9;j,v 00 ,v^)) 



-tt/2 



2 <p(^(0; j, Woo, «*); J, ^oo, 



de. 



Note that in the shallow water case /(f) = 1/v, 



v (~f(v) + Ifv 2 + Xv + n) = (v- v m ) ( |jV + afv, +X)v + 



1 „'2„,2 



' 1 „"2„ 



hence ip is explicitly given by 

<p(v; j, Voo, u*) = |- (t> + f m ) , v m := \ ( \fv* + A + v/ (^j 2 ^* + A) 2 - 2 j 2 /t>* 



First of all, we want to check whether (68)-(71) is an evolutionary system, or in other 
words if the Jacobian matrix M of the mapping 

P := W *) T ^ W(P) := (k, (V), a + j(V), j(V 2 ) + a(V)f 



is nonsingular. This can be done numerically. Next, we can rewrite (68)-(71) in the more 
compact form 

d S W(P) - dy(jW(P) + F(P)) = , 

F(P) := (0, a, A - ja, \a 2 - fif , (j, A, /i) = Voo , u») , 
from which we easily infer the quasilinear form 

(72) M (P)(<9 5 - jdy)P + M X (P)dyP = , 

withMi(P) :=-JacF(P) - j W(P) (1, 0, 0, 0), i.e. 

( -k 



/ 



Mi(j, a, Uoo,i>*) 



-A; 

-j(^)+2jf c 



-A; 



9A 






dX 



Mv)+j(v 2 )) + 



dj 



-a 



dvoo 



\ 







J 






j 2 + p'(Voo) 










\ -(a(V) + j(V 2 }) + 2jVooV* - jvl -a (j 2 + j/(voo))«, p(t>oo) - + j 2 (^oo - ««) / 

By 'change of frame' Y \-> Y + jS, the hyperbolicity of (72) at P = P is equivalent to 
the hyperbolicity of 

M (P)d 5 P + Mi(P)0yP = 

at P = P. Once Mo(P) is known to be nonsingular, a sufficient condition for hyperbolicity 
is that the eigenvalues of Mo(P) _1 Mi(P) have four distinct real parts, because then they 
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must be real and distinct. We comment below on a series of numerical results obtained 
with the shallow water pressure law p(v) = 1/v 2 , in which we have computed mean 
values by using the trapezoidal rule with 10 000 points of discretization, and the Jacobian 
matrices Mo(P) by means of a finite difference method and a discretization step h = 10 -6 . 
In each picture, we have plotted the real part of the eigenvalues of Mo(P) _1 Mi(P) as 
a function of the period S: the modulated equations are hyperbolic for a given k if we 
find four distinct real parts. When two curves collide, there is a set of complex conjugate 
eigenvalue and the system is not hyperbolic any more. 

4.1.2 Numerical results 

We have checked the hyperbolicity of the modulated equations in the cases j = 1 and 
various values of G (0.55,1.2) and in the case j = 4, !)„ G (0.15,0.49). Since the 
hyperbolicity does not depend on a, we have set a = in M and Mi. In both cases we 
found a similar scenario. Let us describe the case j = 1. If > ~ 0.86, the Whitham 
equations are hyperbolic for all periodic waves (see figure [2]). 




Figure 2: Re(Aj), i = 1,2,3,4 as functions of the period H. On the left: j = 1 and 
i>oo = 0.93. On the right: j = 1 and Voo = 0.9. 

If < vlo, one finds a range of periodic wave periods S G (S m ; Em) for which the 
Whitham equations are not hyperbolic, and thus periodic waves are spectrally unstable; 
see figure [3] (left). Moreover, if J$ 0.58, there is one eigenvalue that diverges to ±oo 
which means that there exists S c such that det(M ) = and the modulation equations 
are not of evolution type. If 5 ^ H c , the scenario is similar to the previous case: there is 
a range of unstable periodic waves; see figure [3] (right). 
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Figure 3: Re(Aj), i = 1,2,3,4 as functions of the period S. On the left: j = 1 and 
= 0.84. On the right: j = 1 and = 0.55. 

4.2 Van der Waals type pressure laws 
4.2.1 Parametrization of periodic waves 

We consider in this section pressure laws that are possibly nonconvex, and even nonmono- 
tone. A typical example is the van der Waals pressure law 

RT a 

P{v]T) = - — , 

v — V 2 

which exhibits various types of behaviors depending on the temperature T (compared to 
a/ (bR), where R is the perfect gas constant, and a, b, are parameters of the specific fluid). 
More precisely, defining 

_ 81a _ 8a 
°' 256W?' c ' 27WT 

we easily see that 

1. for T > To, v i — y p(v; T) is monotonically decaying and convex, 

2. for T c < T < T , v i— > p(v;T) is monotonically decaying and admits two inflection 
points, 

3. for T < T c , v i—)- p(v; T) admits one local minimum, one local maximum, and two 
inflection points. 

After nondimensionalization this pressure law reduces to 

7 1 



p(v) 



V — 1 V 



2 
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with 7 := ^-2. As seen above, transition values of 7 are 70 = and 7 r = — . 

' 256 ' 27 

In what follows, we have chosen to deal with the case of a non monotone pressure and 
T = 600K with a and b roughly corresponding to water, so that 7 < 0.275 < j c . 

If we choose such that j 2 < — p^v^), and take as in the previous subsection 
A = — j 2 foo — p{voo), the phase portrait of the travelling wave ODE 

<V)V« + k K '( V ) V < + P( V ) + ^ V + A = °> 

is certainly independent of the function k, but heavily depends on the values of j and v,x, 
when the pressure p is non monotone. 



p 




Figure 4: Left: two fish phase portrait. Right: eyes-and-guitar phase portrait. The 
shaded regions correspond to the domain of hyperbolicity of Whitham's modulation equa- 
tions for the cases we have tested 
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Depending on areas delimited by the pressure curve and by the so-called Rayleigh line, 
of equation y = —j 2 {y — Voo) in the (v,p) plane, two typical phase portraits arise, which 
may be described as follows (see figure Q. 

Two fish There are two disconnected homoclinic loops (the 'fish'), one ending at 

and one with another endpoint, and all the trajectories outside these loops are 
unbounded; thus there are two types of periodic orbits, surrounded by either one of 
the homoclinic ones; 

Eyes-and-guitar There are two homoclinic loops (the 'eyes') ending at the same point 
(maybe t>oo), and a third homoclinic loop (the 'guitar') surrounds them; there are 
three types of periodic orbits, those inside the eyes, and those in between the eyes 
and the guitar. 

A series of numerical investigations in various cases is reported below. 

4.2.2 Numerical results in the eyes-and-guitar case 

We have first considered the case j = 0.023732, = 6.598196 and checked the hyper- 
bolicity of Whitham's equations with a simple form of the capillary coefficient: k(v) = 1. 
Here, we have an "eyes and guitar" type phase portrait. There are three families of peri- 
odic waves. The first one that is 'above' the doubly homoclinic orbit. In this case, there 
exists S c « 559 such that the modulation system is not of evolution type. Furthermore, 
there exists E M w 555.4, such that if 5 > E M and H ^ H c , the modulated equations 
are hyperbolic and thus periodic waves are stable under large wavelength perturbations 
(see figure [5]). On the other hand, we found that the periodic waves we could compute in 
the loops of the doubly homoclinic orbit are stable under large scale perturbations (here 
Whitham's equations are hyperbolic); see figure [6j In the smallest loop, we stopped the 
computations at H ~ 105: below this value, the amplitude of periodic waves is too small 
in comparison to the precision we fixed. Anyway, as explained at the beginning, we expect 
that these periodic waves are unstable for sufficiently small amplitudes, since the Euler 
equations are not hyperbolic at the center point (pressure is nondecreasing at that point). 
We have tested j = 0.032 and the situation is the same. 

4.2.3 Numerical results in the two-fish case 

Here, we have chosen j = 0.0258 and = 1.90285: in this case there are two separated 
homoclinic orbits. The first one is homoclinic to Voo and the second one is homoclinic to 
Woo = 7.57197. In this latter homoclinic, there is a center at wq = 32.49447 from which 
bifurcates a family of periodic orbits. We have checked numerically the hyperbolicity of 
Whitham's equations associated with this family: the modulation system is always hy- 
perbolic (see figure [7|) . 

In the homoclinic orbit to Voo = 1.90285, there is a center at Vq = 3.2089 from which 
bifurcates the second family of periodic orbits. The picture is a bit more complicated: 
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Figure 5: Re(Aj), i = 1,2,3,4 as functions of the period S. On the left: j = 0.023732 
and Voo = 6.598196. On the right: zoom of the previous picture 
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Figure 6: Re(Aj), i = 1,2,3,4 as functions of the period H for periodic waves for j = 
0.023732 and = 6.598196. On the left: periodic waves in the largest loop On the 
right: periodic waves in the smallest loop up to 5 ps 105 (the fourth eigenvalue is real and 
smaller than —10). 



if 90 < H < H m 95.15, the modulation system is hyperbolic. Again, we stopped 
computations at H = 90 in the smallest loop since amplitude of the waves are too small 
(for the precision of computations fixed here): it is expected that for sufficiently small 
amplitudes, periodic waves are unstable since the Euler equations are not hyperbolic there. 
Then if H G (H m ; Em ~ 125.5), the system is not hyperbolic and the associated periodic 
waves are unstable. Finally, if H > E M , the Whitham equations are hyperbolic again (see 
figure |8|). 
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Figure 7: Re(Aj), i = 1, 2, 3, 4 as functions of the period H for periodic waves in the loop 
of the homoclinic orbits to Wqo = 7.57197. Here j = 0.0258 and = 1.90285 




-0.15- 



Figure 8: On the left: Re(Aj), i = 1, 2, 3, 4 as functions of the period S for periodic waves 
in the loop of the homoclinic to v^. Here j = 0.0258 and = 1.90285. On the right: 
zoom of the previous picture on the three eigenvalues with the smallest real part. The 
fourth one is always real 

4.2.4 Conclusion from numerical investigations 

We have found, in the cases where a doubly homoclinic orbit occurs, that periodic waves 
that are inside a homoclinic loop are always spectrally stable under long wave length per- 
turbations or, more precisely, that the modulation equations are hyperbolic. There is an 
additional family of periodic waves for which we observed that there is a critical period 
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S c where the modulated equations are not of evolution type. Moreover, we found that 
there is Em such that Whitham's equations are hyperbolic only if E > 5^. 

If there are two separated homoclinic orbits, we found that periodic waves which 
correspond to the larger solutions pass our stability test whereas there are range of periods 
where Whitham's equations are not hyperbolic for the smaller solutions. 



A A concrete computation 

We derive here the averaged equations associated with Benjamin's impulses pu and vu 
for ([I]) and ^ respectively. This computation is given for concreteness, even though it 
is contained in the abstract computation made in Section [2j Taking the inner product 
of plj) and (42) with (u , Po) and (wo,v ) respectively, then averaging and integrating by 



parts in 9, we receive 

-K (d e (p (u -o))ui + de{\{u$ - a) 2 ) pi + A (d e po) pi) 
+ d T (pou ) + dx(poul) + (p (d x go + K d e B )) = 0, 

k{d 6 (w - jv )w 1 - (jd e w )v 1 - a (d e v )v 1 ) 

+ d s (v w } - d Y (\wl) + (v (d Y p + kd e b )} = 0, 

All terms with index one here above cancel out because of the profile equations (39) and 
(40), which imply indeed that 

d e (po(u - a)) = , d e (l(u - a) 2 ) + A (d e p ) = , 

de(w -jv ) = 0, jd 9 w + a (d e v ) = 0. 

This is straightforward for the equations on the left, and for the ones on the right we 
observe that the second order differential operators A and ao have been defined in such 
a way that 

A (d e po) = dgg , a (dgv ) = %> • 



Finally, we recover the equations in (50) and (|52j) obtained by averaging the impulses' 
conservation laws by checking that 

(po(d x go + Kd 9 B )) = d x (po9o + K (d 9 p ) —(p , Kd e p ) - ^(p ,Kdgp ) 

I de 
(v (d Y po + kdgbo)) = dy^vopo + £(v ,kd e v ) - k(d e v ) —(v ,kd e v ) 

Let us check the first equality, the second one being identical through the symmetry 
(X, K, p , g , B ) -B- (y, k,v ,£, -p , -b ). We have 

dxipogo + K(dgpo)—(p ,Kdepo) - S{po,Kd e po)J - (po(d x go + K d e B )) = 
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/(dxpo)go + d x (k (d e po) J^(Po, Kd g p ) - &(po, Kd e p ) J ^ - K(p d e B ). 
Recalling the definition of go and integrating once by parts (in 9 of course), we get 

/ dS d<§ \ 

({9xPo)go) = ((dxpo) ~^(po, Kd e p ) + (dedxpo) -^-(po, Kd e p )) = 

I d^ \ 

d x (<?(po,Kd eP o)} - ((d x K)(d 6 po)^(po,Kd e po)j = 

I d$ \ I d<§ \ 

dx (&(po, Kdepo) - K (d e po) ^~(Po, Kd e po)j + K d x ((d e p ) -^-(po, Kdepo) j . 

So it just remains to show that 

/ d$ \ 

(p d e Bo) = d x (^(dgpo) —(po,Kd e po)J . 

Once again, this follows from integrations by parts. Indeed, 

/ d 2 <§ ( d<§ 

-((d e po)B ) = (-(depo) q-q~ (po,Kd e p ) (d x po) + (depo) d x \w^-(po, Kdepo) 

( d 2 S \ 
-K(d 2 p ) I —(po, Kdepo) J (dxpo) 

(dS \ ( dS 

— (po,KdePo)J (d x po) + (dePo) 9 X ( ^— (po, Kd e p ) 



J x 



dx ((depo) jr-(po, Kdepo)) 



d<S 

dp x i 

B A convenient structural assumption 

Our purpose here is to check that, under a reasonable structure assumption on the Hamil- 
tonian Jffi the operator srf v behaves properly, and if the kernel of g/^ has the expected 
size, our parametrization hypothesis for periodic profiles is met. 

Structure of the Hamiltonian To go further into the analysis of the abstract equation 
([5]), we need to be more specific about the form of the Hamiltonian Inspired by our 
examples @@, @Q> and @> we write the U-space as R N = R n x 1^^" for some 
integer n, < n < N, require that 

U= (u)' ^(U) = ^(v,u) + ^(v,v x ), 
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and assume that <ffl + c£2 is uniformly strongly convex in both v x and u on the range 
of (U, v x )-values and speeds c under consideration. Note that a simple way to make this 
assumption independent of c is to assume that J -1 has a block structure of the form 

J 1 = 







(N-n)x(N-n) 



as is the case for the Euler-Korteweg system. 



B.l Compactness of resolvents 



We briefly sketch here a proof of the fact that our structural assumption on ensures 
that =2/ (0) with domain # 3 (R/Z; lR n ) x i^R/Z; R( iV ~ n )) have a nonempty resolvent set 
and compact resolvents. In turn, this implies that, for all z/, stf v is a relatively compact 
perturbation of 

Our structural assumptions readily yield 

= jjfea,A<°> a<°> 



S 3 + H 2 kd e - kde^* 2 ~ kd e ^ikd e 



with Si 
S 1 = 



J 3- 





^nx(N-n) 


0(7V-ri)xn 


0(N-n)x{N-n) 



* 


Onx(AT-n) 


0(7V-n)xn 


0(V-n)x(V-n) 



0-3 

with 0i and 03 being uniformly positive definite. Then, standard energy estimates enable 
us to show that 

|Re(z)| ||U|| 2 



(kd e 



C kdeu 



(U, (z - ^(°))U) 
U 



> 
> 



Cllv 



|2 

I #2 



|((/c<9 e ) 3 v,o-i(£;<%) 3 v) + f (£;<9 e u, a 3 /i;<%u) 
-C|Im(z)|[||v||^ /2 + ||u||^ 1/2 ]-C||U|| 2 

where C is a positive constant independent of z. From this we obtain that there exist 
7] > and C > such that if |Re(z)| > 77 [1 + |Im(z)| 6 ] then 

(73) ||U|| M i < C\\(z - s/W)U\\. 

This already shows that, for such a z, (z — g/^) has a closed range and is one-to-one with 
a continuous inverse. To check that the previous range is dense, we only need to examine 
whether the formal adjoint is indeed one-to-one (on smooth functions). This amounts to 
show that 

(74) (-z-A^Jkd )V = 



has no nontrivial smooth solution V. Applying the operator Jkdo to (74), we deduce 
from (73) applied to z and U = JkdgV that kdgV = 0, which in turn implies V = 0, 
because of (74) and the fact that z is nonzero. This proves that for the above 77 > 0, if 
|Re(z)| > 77 [1 + |Im(z)| 6 ] then z lies in the resolvent set of £f(°\ 

Since we have not used any Poincare inequality in our previous arguments, they apply 
mutatis mutandis to srf with domain H 3 (K; R n ) x H 1 (K; M^ - ™)), and show that it has a 
nonempty resolvent set. 
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B.2 Parametrization of periodic orbits 

Let us prove that, under our structural assumption on Jff , there is no restriction in assum- 
ing, as done in Theorem [TJ that Whitham's parametrization by (k, M, P) is admissible. 
More precisely, we are going to show that, if nearby periodic traveling wave profiles form 
a JV + 2 dimensional manifold, and if the generalized kernel of is of dimension N + 2, 
then those nearby profiles are parametrized by (k, M, P). This extends to our Hamiltonian 
framework a result previously shown for parabolic conservation laws by Serre [T9] . 

To set things on a more formal ground, we define on some open neighborhood U of 
wave values (5, c, v(0), v^O), A) the map 

n- u ~^ R2n i 

(H,c,vo,v ,£., A) i — > ([v]f, [v x ]f) 

where [-]q denotes the jump [/]q = f(E) — /(0), and U is the solution of 

E(^ + c^)[U] = A, v(0) = v , v x (0) = v 0>:c . 

We identify in the usual way nearby periodic traveling wave profiles with elements of the 
zero set of 7Z. 

Proposition 2. Assume that U is non trivial and that 1Z has constant rank 2n — 1. Then 
the generalized kernel of is of dimension N + 2 if and only if, up to translation, 
nearby periodic traveling wave profiles may be regularly parametrized by (k, M, P). 

Proof. Our proof is based upon the fact that the dimension of the generalized kernel of 
&tf(°> is the algebraic multiplicity of zero as a root of some Evans function D( ■ ) (see |11|). 
Indeed, viewing spectral problem 

(75) zV = &fV 

for (z, V) = (z, (v, u) T ) as a system of coupled differential equations of third-order in 
v and first-order in u, we may introduce its fundamental solution R(z; •) normalized by 
R(z;0) = Id R 3n xK (jv-n) and define 

D(z) = det([R(z; ■)]!)• 

Then the condition on the dimension of the generalized kernel of reads 

D(z) = az N+2 + 0(z N+3 ) 

for some nonzero a [llj. We want to convert this into some information about profiles 
parametrization. 



Let us denote by V J (z; •) the solution to (75) corresponding to the j-th column of the 
matrix R(z; •), that is V J (,2; •) solves (75) and (y^z; 0), vj(z; 0), v 3 xx (z; 0), u^(z; 0)) T is the 



j-th vector of the canonical basis of lR 3n x MS N n -* . The Evans function is then written 



D(z) 









j v V+2n 






vi] •• 








v 1 1 


. r v JV+2n 


XXI 


[ XX 


u 1 ! ■■ 


. \ u N+2n 
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where we have dropped the marks and H on jumps. To go further, using our structural 
assumptions, we writ^] 



3d x A, A = £ 3 + £ 2 <9, - dJZ - 



with Si = £*, S 3 = £*, 
S 1 = 





Onx(iV-n) 


0(7V-n)xn 


O(JV-n) x(iV-n) 



* 


Onx(iV-n) 


0(7V-n)xn 


O(JV-n) x(iV-n) 



0-3 



o"i and CT3 being uniformly positive definite. Integrating (75) from to 5 yields 



z 



*i(0)[v£J+ *[vj] + *[v^] + *[u 



a 3 (0)[u^]+ *[v^ 

Therefore, up to a nonzero multiplicative constant, D{z) is also written 



[v 1 ] 
[vi] 



N+2nl 



^V 1 



Now, corresponding to the impulse equation, we also have 

z lr (v « } = a * ( n " Hess( ^ + ^ )(VJ) + vj ' E( ^ + - d{Jf au ~^ (v ° } 

= d x (U- AW + EjJJ. • Vj + ((S 2 - S*)U X - SiU.J • V^) 

where the convention is as before that linearization and derivatives are taken at U. By 
integrating the relation here above, we obtain 



-as 

dU a 



(V4) = (aiv,)(0) • K] + * M + * z / V. 

Jo 



We still need to check that it is not a trivial relation. But, since U is non trivial, there is 
a point where is nonzero, otherwise v would be constant thus u and U would also be 
constant, a contradiction. Then, since assumptions of the proposition and terms of the 
equivalence we are currently proving are invariant by translation, we may assume that 
v^O) is nonzero. Now let us pick £ such that the £-th component of cr 1 (0)v a ,(0) is nonzero 
and, for any V = (v, u) T e R N = R n x R( N ~ n \ denote by v* the vector of R^'^ obtained 



3 We warn the reader that, because of k factors, these notations are not compatible with the ones of 
the previous section of the present Appendix. 
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from v by deleting the £-th component. Then, up to a nonzero multiplicative constant, 
D(z) is 



y iV+l 



[ (vj) J 
rl as /yi- 

JO dU„ \ v a. 



r N+2n] 



N+2n\ 



10 8U a 

Jo-V 1 



[(v 



C=l -yN+2n 

10 ' Jo 

Up to a change of basis we may assume that V 1 (0; ■) = U T , V J '(0; 



U; 



2n + 1 < j < iV + 2n and (V^O: 



>i-2n 



for 



V (0;-)) is a basis of the linear span of 



{U 

D(z) is written up to a nonzero multiplicative constant 



(v )l 



U( vo v }. With this choice, after setting V 1 = V z (0; 



M+2 



v 1 ] 


[v 2 ] 




(vi)J 






Jo 8U a \ V 


Jo 8U a \ v aJ 


f=L 83 (TjN+2n 
Jo 8U a \ V a 


Jo-V 1 


Jo-V 2 





+ 0{z 



N+3\ 



Since V 1 satisfies XJ_ X = it differs from U c by an element of the kernel of srf 

which is spanned by (V^O ;•),-•• , V JV+2n (0 ; •))• This yields that, up to a multiplicative 
nonzero constant and an additive remainder 0(z N+3 ), D(z) reads 



z 



N+2 



or 



v N+2 



vjo) 

(vj,,(0) 
1 

<2(U)(0) 
U(0) 



[v c ] 

[((V*)c)x] 

/o-U c 



[v c ] 

K(v*)c)*] 
j] 

Jo-g((U c 



[v 2 
[(v 

Jo dU a v Q - 
Jo-V 2 



[V 2 (0;-)] 
[(v 2 ),(0-)] 


J -^(V 2 (0;-)) 



Jo" V2 (0-) 



[v^ 2 ™] 

f ~ 83 (yN+2r 
Jo (9f7 Q v " 



[v* +2 "(0;-)] 
[(vf +2 "),(0; 




•)] 



ff(vr(o;-)) 

J?V"+*»(0;.) 



We are ready to complete the proof by observing the latter determinant. Indeed 
our assumption on TZ implies that the (2n — l)-st rows of the above matrix are linearly 
independent. Furthermore, the kernel of the corresponding linear map is the tangent space 
at U of the profiles manifold (profiles being identified when equal up to translation). Thus 
the differential map of U (H, =2(U), f~ U) is invertible on this tangent space if and 
only if the above determinant is non zero. Consequently, this map is full-rank if and only 
if the generalized kernel of g/^ is of dimension N + 2. □ 

Note that, by introducing Floquet exponents in previous arguments as in analogous 
computations in [19] , one may obtain an alternative proof of Theorem [TJ 



48 



C Galilean invariance 

In order to check that the hyperbolicity of the Euler-Korteweg modulation equations 
(53) does not depend on a, it is convenient to rewrite those equations in a form that is 
similar to (66) for (55). Substituting ga + j for the mean momentum (po^o) m (53), and 
manipulating the remaining mean values as in the proof of Theorem [2j we receive the 
system 



(76) 



' d T K + d x (aK) = 0, 
d T g + d x (Q(T + j) = 0. 



<9- 



T 



ga + j + D 
Q 



+ <9 



x 



l(ga + j)(ga + j + 2D) 



d T (ga + j) + d x 



2 e 
(ga + j) 2 + 2jD 



+ 



Q 



+ £g + KB — E 



J 



together with the generalized Gibbs relation (62) dE = gd£> + QdK + - dD. Then, we 



easily check that ( 76 ) is invariant by the Galilean transformation 

(T, X, K, g, a, D) (->■ (T,X — aT, K,g,a- a, D) 

for any a. Since this transformation leaves invariant all the 'thermodynamic' variables 
(K, g, D, g, 9, j, E), it leaves (76) invariant just because it does so for the reduced system 

' d T K + d x (aK) = 0, 
d T g + d x (g(r) = 0, 

d T a + d x -cr H cr 

V. 2 Q 

k Mo*) + d x (Qo- 2 + 2ja) = 0. 
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